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Abstract 

The renormalized charge of a simple two-dimensional model of colloidal 
suspension was determined by solving the hypernetted chain approxima- 
tion and Ornstein-Zernike equations. At the infinite dilution limit, the 
asymptotic behavior of the correlations functions is used to define the 
effective interactions between the components of the system and these 
effective interactions were compared to those derived from the Poisson- 
Boltzmann theory. The results we obtained show that, in contrast to the 
mean-field theory, the renormalized charge does not saturate, but exhibits 
a maximum value and then decays monotonically as the bare charge in- 
creases. The results also suggest that beyond the counterion layer near to 
the macroion surface, the ionic cloud is not a diffuse layer which can be 
handled by means of the linearized theory, as the two-state model claims, 
but a more complex structure is settled by the correlations between mi- 
croions. 

PACS: 82.70.Dd, 61.20.Gy, 61.20.Qg 

1 Introduction 



The omnipresence of colloidal suspensions in chemistry and biology, and the 
variety and tunability of their particle interactions explains their importance 
both in industrial applications and fundamental research. A complete statistical 
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description of this kind of systems is an extremely complex task inasmuch as that 
description must include the interaction between all of the species constituting 
the system: colloids, suspended ions and molecules of the solvent. However, as 
a result of the high asymmetry of mass and charge, one hopes that a suitable 
approach is obtained by considering the system, not as a mixture, but as a 
"new" monodispersed system which involves only colloidal particles (macroions) 
interacting through an effective potential [U [2] . 

In the framework of the approaches based on linearized Poisson-Boltzmann theo- 
ry, it is found that the effective interaction between two macro-ions exhibits a 
Debye-Hiickel-like (DH) behavior due to the screening effect produced by ions 
present in the surrounding medium. As the charge of macroions (Q) increases, 
the linear approximation becomes inadequate and it is necessary to use other 
approaches. Nevertheless, these approaches show that at large distances the 
effective interaction continues to have a DH-likc behavior but it is necessary to 
replace the charge Q by an effective or rcnormalizcd charge Qren- Then, the 
concept of charge renormalization simplifies the description of colloidal inter- 
actions and gives us some insights about the phase behavior of the colloidal 
suspensions. However, it remains problematic to find methods that allow their 
evaluation for arbitrary conditions of charge and concentration. 

In this work we evaluate the renormalized charge for a two-dimensional (2D) 
Coulomb system. For this purpose, we solve numerically the Ornstein-Zcrnike 
(OZ) equation in the scope of the hypernetted chain (HNC) approximation 
and we use the calculated potential of mean force as the effective interaction 
between the components of the system. The renormalized charge, numerically 
evaluated, is then compared to the one predicted by theoretical models based 
on the Poisson-Boltzmann (PB) theory. The results obtained show that, in 
contrast to the mean-field theory, the renormalized charge does not saturate, 
but exhibits a maximum value and then decays monotonically as the bare charge 
increases. The results also suggest that beyond the counterion layer near to the 
macroion surface, the ionic cloud is not a diffuse layer which can be handled 
by means of the linearized theory, as the two-state model claims, but a more 
complex structure is settled by the correlations between microions. 

We use a two-dimensional system due to two main reasons. In the first place, 
classical two-dimensional Coulomb systems with logarithmically interacting par- 
ticles are simplified models which keep the Coulomb nature of the particle in- 
teraction. Besides, by using field theoretical tools, some of these models admit 
an exact solution, which makes it possible to judge the quality of approximate 
formulations. On the other hand. 2D Coulomb systems can be used to model 
stiff rod-like polyelectrolytes, which interact with the ions present in the solvent 
via a logarithmic potential. 

This article is organized as follows. In section 2 we describe briefly the system 
considered. Section 3 is devoted to explain the numerical method used to solve 
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the OZ and HNC equations. The central part of this work is presented in 
section 4, where the numerical results of the HNC calculations are displayed 
and discussed for two models of 2D Coulomb systems: the guest particle system 
and a simple model of colloidal suspensions at the infinite dilution limit. Finally, 
some conclusions are summarized in the last section. 



2 Guest charge in an electrolyte 

In order to take the first steps toward a solution of the renormalized charge 
problem in colloidal suspensions at the infinite dilution limit, we consider the 
properties of a "guest" charge immersed in an electrolyte, which is modeled as a 
two-dimensional two-component plasma (2D-TCP). We consider here a system 
composed by a single central guest particle with charge Q = Ze and radius tro 
immersed in a 2D-TCP formed by negatively and positively hard-disks ions of 
radius a± ~ a and charges q± = ±z±e (with z± > 0). The bulk densities of 
positive and negative ions are denoted by n_|_ and n_, and the total density 
is 71 = ?i+ + 7i_. The electrolyte is globally neutral, q+n+ + q-u^ ~ 0. The 
interaction potential between disks i and j, separated by a distance , is given 



where /3, z^, and T = /?e^ are, respectively, the inverse temperature, the valence 
of the charge i and the coupling (strength) parameter. L is an irrelevant length 
scale used to fix the zero of the Coulomb potential. At the limit of point-like 
charges, CTq = and cr = 0, the systems is stable against the collapse of positive 
and negative particles provided —2/2+ < ZT < 2/2_ and z+z_r < 2. 

From field theoretical arguments in the limit CTo — > and a = 0, the asymptotic 
large-distance behavior of the effective interaction between the guest particle and 
the TCP charges is related to the large-distance behavior of the corresponding 
two-points correlation function of the exponential fields associated to the 2D 
sine-Gordon theory [3l[4]. According to this correspondence, the asymptotic 
behavior of the effective interaction can be exactly evaluated and this result indi- 
cates that the electrostatic potential produced by the guest charge has the same 
asymptotic behavior as the one evaluated through the DH approach provided 
that both the inverse decay length and the charge are adequately renormalized. 

On the other hand, the short-distance asymptotics of the density profiles of the 
colons and counterions near the guest charge has been studied for the symmetric 
case, i.e. z± = 1 [5l[6]. In this case, it has been shown that the colon density 
profile exhibit a change of behavior if the guest charge becomes large enough 
{ZT > 2— r). This behavior has been interpreted as a first step of the counterion 
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condensation (for large coulombic coupling); the second step taking place at the 
usual Manning-Oosawa threshold {ZT = 2). 

Naturally, the latter effect is related to the settlement of the renormalized 
charge. In this way, it can be stated that this model encloses the most im- 
portant facts involved in the non-linear screening phenomenon. In order to 
compare these exact results with the numerical ones in more general conditions, 
we have calculated the correlation functions between the different components 
of the system by solving the OZ and HNC equations as described below. 



3 Numerical method 



As usual, we denote by gij{r), hij{r) = gij{r) — 1 and Cij{r) the radial distribu- 
tion function, the total correlation function and the direct correlation function, 
respectively, between two components i and j of the system separated by a 
distance r. The subscripts i and j can be +, — , or 0, denoting respectively 
a positive ion, a negative ion, or the guest charge. The OZ equations for a 
multicomponcnt system reads 



hij 



ki - i"2|) = Cy(|ri - r2|) + ^rife / /iifc(|ri - rsDcfcj (|r3 - r2|) dra (2) 

7. 



where Uk the density of the species k. In the infinite dilution limit considered 
here, no = 0, the OZ equations describing the correlation functions between ions 
of the TCP do not involve the guest charge. By defining 7y (r) = hij (r) — aj (r) , 
these equations can be written in the Fourier space as 

„ (1 - Djk)) c++{k) + n_ (c^_(fc) - c++(fc)c__(fc)) 

^ m 

„ (1 - D{k)) c^^jk) + n+ {cl_ik) - c++(fc)c__(fc)) 

= m 

where the hat denotes the Fourier transform of a function, and D{k) = (1 — 

n+c++(A:))(l — n_c (k)) — n+ri-c^ (k). The guest charge-ion correlation 

functions, in the infinite dilution limit, satisfy 

„ co+{k){l-n-c — (k) - D{k)) + n^co-{k)c+-{k) 

70+ (fc) = (6) 

, co-{k){l-n+c++{k)~D{k)) + n+co+{k)c+^{k) 

10- {k) = . (7) 

In order to solve equations ©-([T]) it is necessary to split functions 7(r) and 
c(r) in short- and large-distance components. Since h{r) y c(r) —j3v{r) 
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when r — > oo, it may be written 

i;(r) = i;(-^)(r) +w(')(r), c{r) = c^'^r) - I3v^^\r), -/{r) = j^'\r) + (3v^^\r) 

where the large-distance component of the potential is chosen to be v^^'>{r) = 
(1 — e~°'^)v{r) with a a cutoff parameter. In this way, the corresponding HNC 
closure relations are defined through 



r < a. 



^exp(7l;^(r)-/3T;J)(r)) 



r > (T, + cr ," 



c. 



^'\r)=g,,{r)-l-^lf{r) (8) 



with v^^f (r) the short-range part of the Coulomb potential. 



Equations ©-([S]) were solved by using a direct iterative method. In order to 
improve convergence, the new input function for iteration i -|- 1, c|^~'"^^(r), was 
evaluated by using a linear relaxation (Broyles' mixing of the old input function 
c[l^ (r) and the output function CQ^\(r)) 

^'^(r) = ci^ir) + f^(^,(r) - c^ir)) (9) 

where (0 < ^ < 1) is the relaxation parameter. The numerical evaluation 
of the two-dimensional Fourier transform was carried out through the quasi- 
discrete Hankel transform method [7| , by using = 6000 grid points and max- 
imum radius 10 < kR < 16, with k the inverse Debye length (k^ = 
The iterations were stopped when the larger difference between successive c(r)'s 



T = max|citii)(r)-c«M|<10-^ (10) 

Finally, the numerical accuracy was monitored by checking both the elcctroneu- 
trality condition and the Stillinger-Lovett screening sum rule [5] . In all the con- 
sidered cases, the numeric results for these conditions differ at most by 0.01% 
with regard to the theoretical values. 



4 Results 



In order to calculate the renormalized charge Zren we follow a similar procedure 
to the described one in [Hj . Let us recall that under the DH approximation the 
potential of mean force w^^{r) between the guest charge and an ion of type 
s = ± is given by 

-/3<"M--^g-TA'o(-) (11) 
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with Kq and Ki the modified Bessel functions of order and 1, and dc = o"o + (7. 



Once the correlation functions go±(r) are iteratively evaluated from the OZ 
equations, the potential of mean force wo±ir) is obtained from the definition 
go±i'r) = e~^'^°^^'^\ We assume that, at large distances nr ^ 1, the poten- 
tial of mean form has a DH-likc functional form, i.e. w{r) cx Ko{Krenf) ~ 
^/tt exp{— Krenr) / \/ {2Krenr). Bcsidcs, we can also assume that all correlation 
functions fall off in the same way at large distances. These assumptions are 
supported in the exact solution for point-like particles, which indicates that in 
the stability region the correlation functions exhibit a DH behavior. However, 
it is necessary to renormalize not only the charge Zren but also the decay length 

Kren [31 ■ 

In this way, we may suppose that 

- Pwosir) = \ngos{r) - ^o.^^^^^^^ (12) 



where 



Aos = I TZreiiZs ^^^^ 

Ki{KrenCrc) 

Then, the parameters Aq^ and Hren can be obtained from a linear fit of the 
functions 

Los{r) =ln{^/r\\iigos{r)\) = In {\Aos\) - Krenr (14) 

for r values large enough. We select the range of r to be used in the fit through 
the condition 10~^ < | log f/os (j') | < 10~^. We choose this range because the 
values of the fitted functions are two orders of magnitude greater than the 
tolerance and the numerical uncertainties (see eq. pU]) ). In some cases, the 
correlation functions exhibit an oscillation and then it also becomes necessary 
to impose empirically a minimum value for nr. 



4.1 Guest charge 
4.1.1 Symmetric TCP 

We consider all particles as equally sized hard-disks satisfying nao — na = 
10^^, and the ions have charges z± = 1. In figure [T] are shown the correlation 
functions between the guest charge and the ions of a symmetric TCP for different 
values of Z > and the coupling parameter T. Qualitatively, these profiles 
arc consistent with the common image of an ionic cloud (mainly formed by 
countcrions) surrounding the central charge. From the figure, it is clear that for 
r = 0.2 the correlation functions go±i''') behave monotonically for all displayed 
values of ZT. However, when F = 0.4 and ZT = 3 the correlation function 
30+ (^) exhibit a peak at Kr « 0.1, which rise up as a consequence of the increase 
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of go- jointly with more marked microions correlations, as it can be seen in 
figure [2l 

As stated before, the renormalized parameters were evaluated from a linear fit 
(in a least squares sense) of the functions LQ±{r) similar to the ones shown in 
figure [3l whose non-linear behavior extends roughly from nr = to nr ^ 1. 
From now on, a computed renormalized parameter is equivalent to the mean 
value of the two ones obtained from both Lo+{r) and Lf)-{r) taken indepen- 
dently. Note that the curves were calculated at ZT = 3 which correspond to 
strong coupling regime between the guest charge and the microions. 

We explored a range of couplings O.f < F < 0.5. Within this range, the renor- 
malized inverse screening length K^en is found to be systematically smaller than 
K, and it is a decreasing function of the coupling F, in agreement with the exact 
results [3]. However Kren/n only differs from 1 by at most 2% for the largest 
coupling F = 0.5. 

Figure [4] displays Zren at different coupling parameters. In contrast to the 
PB prescription, Zren is not a monotonically increasing function of Z, but it 
displays a maximal value after which it drops to a point where the implemented 
iterative method fails to converge. This behavior qualitatively resembles the 
one predicted by the exact solution when the microions and the guest charge 
are considered as point-like particles (cr = 0) and a small hard disk (kctq ^ 
0), respectively. According to this result, the renormalized charge exhibit a 
maximum at ZmV = 2 — F/2 and decrease until the collapse point ZcT = 2 at 
which Zren = — 1, as indicated by the dotted line in figure [H This behavior 
reflects the effect of the counterion condensation, i.e. the renormalized charge 
associated with the collapse value Zc of the bare charge is identical to that of 
the bare charge Zc — 1 because of the condensation of one counterion from the 
TCP onto the guest particle [31 [3] . By increasing Z further more countcrions 
become condensed and the previous description predicts a renormalized charge 
oscillating between two extremes. 

As it can be noted from figure IH as F decreases, Zren moves closer to the mean- 
field PB theory predicted values; however, the renormalized charge predicted 
by PB theory is always greater than the one computed from HNC approach. 
Naturally, this behavior can be ascribed to the correlations between the com- 
ponents of the system and in a minor grade, to the finite radius of the ions, as 
it will become clearer below. When the coupling becomes stronger, the HNC 
collapse point is closer to the one predicted for point-like particles, i.e. ZT = 2. 
Although, in the other cases the hard-disk potential inclusion allows us to go 
beyond this point, at F = 0.5 we found that the coupling between the guest par- 
ticle and the countcrions is strong enough to provoke the divergence of the HNC 
method at smaller values of ZT. This effect can be related to the short-range 
structure of the correlations between the guest particle and the microions. 
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By using the sine-Gordon model for the TCP, it has been argued that a change 
in the short-range behavior of the coion density profile can be interpreted as a 
"precursor" of the (Manning) counterion condensation [5]. When ZT increases 
above 2 — F {ZT < 2, F < 2), the coion cloud n+{r) — n+go+{r) shows a change 
of behavior near the guest particle. This fact is reflected in the parameter ao± 
describing the potential of mean force for short-distances defined below. Near 
the guest charge this latter must correspond to the Coulomb potential, i.e. 

- /3wos(f) = ln.gos('') ^ aoslnr (15) 

where aos is a constant depending on the bare charge and temperature. Analo- 
gous arguments, involving the electrostatic potential ip{r) calculated from the 
exact PB equation, establish the change in the short-range behavior of ip{r) as 
a fingerprint of the counterion condensation phenomenon |10j . 

From the field theoretical approach [SJ [S] , it is expected that 

zr<2-F 

"0+ = <! „ ^ zr > 2 - F " " ^ ' 

i.e. at Zr > 2 — F the colons interact with the guest charge as if the latter one 
carries a charge Z = 2/F — l = Zc— 1 independent of Z, whereas the countcrions 
always interact with the bare charge Z |5] . Figure [5] shows the values of aos 
evaluated from the short-range values of gosif)- In order to evaluate them, we 
carried out linear fits of higQs{r) according to (jl5p . The comparison between 
the numerical results at F = 0.5 to the ones given by p6)) indicates that the 
short-range behavior of go+(r) when kct = 10^^ is qualitatively consistent with 
that of a system of point-like charges. However, from the figure it can be noted 
that as zr — *■ 2 a change in go-{r) is also found, as it was suggested in [5]. This 
fact indicates that as the coupling increases, the counterions "see" the guest 
charge with a smaller charge than Z. At this point, the effect produced by the 
correlation between the counterions jointly to the hard-disk potential provoke 
the increase of distance between the guest particle and the countcrions. Also, 
as a consequence of this change at large enough coupling, the maximum value 
of Zren do not rcach the saturation value predicted by PB theory, which is 
indicated in figure |4] through the dashed line at ZT > 2. 

On the other hand, the charge accumulation around Z can be also measured 
through the integrated charge distribution 

27r 

where p{r) = (nz-^-zJ) ((7o+(^) ^ .9o-('')) + •^-)- ^(^) indicates the overall 
electric charge found within a disk of radius r and can be used to define the 
cutoff distance, the so-called Manning radius r™, which separates the condensed 
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counterions from the uncondensed ones. According to a geometrical construc- 
tion, the cutoff distance is identified by an inflexion point in P{r) when it is 
displayed as a function of Inr [TOl El [Hj , i.e. 

P(r) 

Figure [6] shows P{r). As expected, at larger ZT, the short distance behavior of 
P(r) is dominated by the counterion density. Besides, as demonstrated before 
in the third panel of the figured] for which ZT = 3 and F = 0.4, the onset of 
a peak in colon profile at high enough ZT and F docs mark the extension of 
condensed counterions layer and give rise to inflection points in the integrated 
charge distribution. 

The locus of points satisfying Pyyljin) = is plotted in figure [71 This figure 
shows that for all values of ZV an inflection point beyond the Debye length 
is found. This point represents a characteristic change of convexity of P{r) 
(as a function of Inr) induced by screening not by the condensed layer [12]; 
the inflection points associated with this layer appear when ZT > 2. As it 
was mentioned above, for point-like particles it is expected that counterion 
condensation begins at 2— F; however, as consequence of finite radius of particles, 
the onset of inflection points for different values of F takes place at larger values, 
as can be inferred from figure [7] It is important to note that smallest values of 
Tin arise as an effect of the hard-disk radius, as becomes clearer from the inset. 

As figure [7] shows, increasing F produces a thicker layer of condensed coun- 
terions, that is, the counterion condensation is enhanced as a consequence of 
increasing microions correlations. If the rcnormalizcd charge is considered as a 
result of the condensation phenomenon, then a smaller value of Zren^ than the 
one predicted by PB theory is expected. It is indeed the case, as the figure [4] 
displays: at high enough ZF, the renormalized charge decreases as the bare 
charge increases because r™ becomes larger. 



4.1.2 Asymmetric TCP 

In order to study the effects of the charge asymmetry on Zren, we carried out 
a similar analysis to the one presented in the previous section when particles 
forming the TCP are characterized by z+ = 1 and z_ = 1/2. If all particles were 
point-like, the stability region would be restricted to —2 < ZT < 4. We choose 
to work with these "valences" because, some exact results exist concerning the 
pair correlation functions of this kind of TCP (TSj [14]. In this section, we 
consider both negative and positive values of Z, so that, in the cases when 
Z > and Z < we may refer to them as 1:^ and i:l cases, respectively. 

Figure [8] shows the correlation functions between the microions at two different 
values of F. Due to our approximation, in the infinite dilution limit the correla- 
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tion functions gss'{r) are independent on Z. As it can be seen, at F = 0.2 and 
r = 0.7 the agreement between the numerically evaluated correlation functions 
gss' (r) and those predicted by field theoretical arguments in [13] (valid at large 
distances) is fairly good at large enough distance (rer > 0.5 and nr > 1.5, re- 
spectively); this is particularly the case for (r) and g (r). Additionally, 

the electroneutrality condition is fulfilled for this TCP. It was found that the 
zero moment defect was less than 2 x 10"'^ for all values of Z and T reported in 
this section. 

Figures [5] and [TO] show plots of the correlation functions between the guest 
particle and the microions at different values of the bare charge for 1:^ and ^:1 
cases, respectively. These curves were calculated for F = 0.2 and F = 0.7. As it 
can be expected, the correlations between the guest particle and the countcrions 
exhibit larger values when the latter ones have higher valence. It should be noted 
that the curves behave monotonically and do not cross for any coupling in case 
l:i; the same happens at F = 0.2 in i:l case. However, for the latter case at 
F = 0.7 it is observed at larger Z, that the counterion profile near the guest 
particle (which is associated with go+(r)) falls off much faster; also, the colons 
profile shows a peak at Kr « 0.1. As in the symmetric case, these facts suggest 
the presence of a counterion shell strongly bounded to the guest particle, which 
considerably screens the repulsion between the guest charge and the colons, and 
consequently plays an important role in the renormalized charge Zren setting. 

As in the symmetric case, the renormalized parameters were evaluated by a 
linear fit of the functions Lo±{r) as those shown in figure [TTl where the curves 
displayed were evaluated at F = 0.7 and for different values of bare charge Z. 
In addition to the condition lO^"' < | log f/os ('') I < 10^^, in some cases it was 
mandatory to introduce empirically a minimum value of r in order to select an 
adequate range to be fitted. Such a situation was found particularly at large 
coupling and for Z < 0, where the correlation function between the guest charge 
and the colons exhibits a pronounced oscillation as it is shown in figure [Tl] for 
Z = -4. 

Within the range of couplings explored, 0.1 < F < 0.9, the renormalized inverse 
screening length Kren is found to be systematically larger than k, and it is an 
increasing function of the coupling F, in agreement with the exact results [5]. 
Notice that this is the opposite behavior than the one found in the symmetric 
electrolyte. The ratio Kren/i^ only differs from 1 by at most 3% for the largest 
coupling F = 0.9. 

In figure [T^ the renormalized charge resulting from the fitting procedure is il- 
lustrated as a function of Z. The global behavior of Zren shares much of the 
features described in |14| . where an analysis based on the integrablc complex 
Bullough-Dodd (cBD) model was used. The most remarkable outcome corre- 
sponds to ^:1 case at F > 0.8. In such a situation, we observe a very strong 
neutralization of Zren as \Z\T increases beyond 2, which gives rise to the charge 
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inversion phenomenon, i.e. the renormahzed and bare charges have opposite 
signs. 



Because of the asymmetry, the correlations between microions provoke an extra 
attraction of the counterions to, and repulsion of the coions from, the guest 
particle. This causes the accumulation of a large number of counterions around 
of the guest particle, as shown in figure 1101 whose charge exceeds the amount 
necessary to neutralize Z . As it was mentioned, charge inversion is predicted by 
an exactly solved 2D model based on cBD theory [T3]. However, it has also been 
reported in systems having other symmetries by using the HNC approximation 
[9l [15] . The decisive characteristic of the systems in which this process takes 
place is the asymmetry (either of size or charge) among coions and counterions. 
Particularly, when the ions are equally sized particles, it seems necessary that 
the counterions have larger valence (in absolute value) than the coions, as in 
the present situation. 

The charge inversion phenomenon involves microions correlations and therefore 
the mean field PB treatment is not able to predict it. In addition to the charge 
inversion, in the case \:\ the "overshooting" phenomenon predicted in (TB], 
i.e. effective charge becomes larger than the bare charge for intermediate values 
of the latter one, is also observed. In such situation, coions are repelled further 
away in the double layer than they are in the symmetric case. The inset of 
figure [12] illustrates the comparison between HNC with PB results; the latter 
ones were obtained from (<to = cr = 0) |10| 



Equation ([T9)) is strictly valid when — 1 < ZY < 2, i.e. if ZT does not exceed 
the Manning thresholds for counterion condensation, otherwise Z^erX saturates 
to and respectively p^ . 

As in the symmetric case, condensation thresholds mark the change of behavior 
of the correlation functions at short distance [Hj. In figures [T^andlTilthe short- 
range correlation functions and the ao± factors are plotted (see equation (|15p ). 
As before, the q;o± factors determine the strength of the potential of mean force 
wo±(f) between the guest particle and the microions at short distance. The 
curves were evaluated at Z = —3.6 and for different values of the coupling 
parameter F. From the figure we can observe that there exists a point from 
which the interaction between the guest particle and the coions (i.e. negative 
ions) becomes less and less repulsive. As before, that change can be interpreted 
as a first step in the counterion condensation [T3|. At F = 0.8 that change is 
much more dramatic to the extent that the interaction begins to be attractive. 
Subsequent to the sign change in the short-range of WQ-{r), the inversion of the 
renormahzed charge (which is evaluated at large distance) takes place. On the 
contrary, the interaction between the guest charge and the counterions is hardly 
modified as \Z\ increases, i.e. at all coupling parameters, the counterions seem 




(19) 
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to "see" just the bare charge at short distance. 

In order to study the thickness of the condensed layer we use again the inflection 
point criterion for the integrated charge distribution P{r). Figure [T51 shows 
the locus of the inflection point for both i:l and l:i cases. In the first place, 
for the case l:i it is found that inflection points are located beyond the Debye 
length when T < 0.8. In this way, for the considered range of Z > the 
renormalized charge is well described by the PB approach at coupling strength 
below 0.8. On the other hand, when F > 0.8 and Z > 4.5 a new inflection point 
is detected at ~ 0.1, i.e. at these bare charge and coupling parameters 
values the correlations give rise to the counterion condensation. Note that at 
such conditions the renormalized charge is close to reach its maximum value, as 
figure [T^l suggests. 

With regard to the ^:1 case, a richer scenario is found. For a large enough 
coupling, two infiection points are found, the smallest of them corresponding to 
rm- As in the other cases, the onset of is associated with coupling conditions 
at which the maximum (absolute) value of Zren is reached. Upon a further 
increase of |Z| these infiection points shift one toward the other and finally 
coalesce and annihilate, as it can be seen in figure 1151 In addition, the bare 
charge at which such annihilation takes place, Z*, is related to the onset of a 
peak in the colons density, i.e. g^^{r). 

Within the cell model of polyelectrolytes, that annihilation occurs when the 
typical salt screening length (Debye length) interferes with the size of the con- 
densed counterion layer and therefore it becomes unnecessary to distinguish 
between condensed and uncondensed counterions |12j . Moreover, it has been 
argued about the merit of the inflection point criterion, which unifies the phe- 
nomenology of an infinite dilution/finite salt system (as the studied one here) 
and that of the a finite density /vanishing salt cell model [inillTj. However, our 
results suggest a different situation: because of the microions correlations, for 
Z > Z* a. second layer mainly constituted by colons starts to be formed at some 
distance of the guest particle. 

On the other hand, according to the two-state model, the distribution of mi- 
croions around the macroion can be divided into a condensed counterion region 
and a free population; the latter one being considered within the linearized PB 
theory. So, the role of a condensed layer is to impose a new boundary condition 
to the electric potential and consequently fixing the value of the renormalized 
charge [T7]. However, since there is no inflection point of P{r) for Z > Z* ^ 
we may not differentiate such populations at large enough coupling (i.e. strong 
microions correlations). In fact, the renormalized charge and subsequently the 
charge inversion are determined by a "layered" ionic cloud, as the insets in fig- 
ures [To] and [15] suggest. An additional discussion about that "layered" cloud 
will be carried out later. 
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In this way, the process behind charge inversion appears to be a similar one 
to the secondary cause of the "Mechanism I" proposed in [15]. Since the net 
charge around the guest particle has opposite sign to that of Z, coions at some 
distance from the guest charge are attracted electrostatically towards the guest 
particle, while the counterions are repelled electrostatically from it, leading to 
the existence of a range of r-valucs at which the mobile ion system is locally 
opposedly charged, as figure [TBI corroborates for T ~ 0.8. 

4.2 Colloidal suspension at infinite dilution limit 

As it was previously mentioned, 2D Coulomb systems can be used as a first 
approach to model some rod-like colloids and polyelectrolytes. Many synthetic 
as well as biological examples can be classified within this category. Depending 
on the number of charges per monomer unit bound to the backbone, the linear 
charge density A can be tuned, so that typical dimensionless values comprise 
Id = \M^B = 3 — 6.6, where Is is the Bjcrrum length. 

In order to study the effective charge in a little more realistic model of a cylindri- 
cal macroion immersed in an electrolytic solution, we follow similar assumptions 
to those given in [9] and [14]. At first, wc consider a long charged cylindrical 
macroion with diameter ag = 20A immersed in an electrolytic solution composed 
by ions with diameter a = 5A, all of them interacting through a logarithmic po- 
tential. Temperature is fixed through the Bjerrum length Ib = 7.2A, which 
corresponds to aqueous suspensions at T = 293 K. In order to study the ef- 
fect of salinity wc consider both symmetric (1:1) and asymmetric (1:2 and 2:1) 
electrolytes at different concentrations n. 

By assuming the DH-like decay of the correlation functions <?o±('') at large 
distances, we have obtained for the symmetric case the renormalized charge 
exhibited in figure [T7J We found a qualitative behavior that resembles the 
one shown by the guest charge system. Again, the most significant difference 
between the HNC and analytical PB results [10] is the presence of a maximum 
value of £,ren followed by a monotonic decay in the HNC outcomes. 

A similar dependence has been described for the zeta potential of a charged 
cylinder [l^ and for the renormalized charge within a spherical cell model of 
ionic condensation [19l [20] . By using the density-functional theory based on 
a weighted (i.e. non-local) density approximation, the monotonic decrease of 
S^ren as ^ incrcases can be explained in terms of the excess chemical potential 
l^excii")- As ^ increase, the number of counterions is also increased due to 
the electroneutrality condition and therefore fiexcif) near the particle surface 
decreases, resulting in an additional number of counterions attracted towards 
the surface |19| . In this way, ^ren must decrease when ^ surpasses the value 
where the maximum ^rem which is as well related to the plateau predicted by 
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PB theory, is reached. 



As before, the extent of the condensed counterion layer can be measured by 
using the integrated charge distribution (fTT)) . In figure [TH] inflection points of 
P(r) are plotted as a function of the bare linear charge density for different salt 
concentrations. At high enough density {n > 0.05M), the counterion conden- 
sation appears to occur at the macroion surface and involves a single layer of 
counterions. This layer is moderately bound to the charged cylinder, as it can 
be concluded by comparing the order of magnitude of go-{ac) at n = O.OOIM 
to that at n = O.IM. Beyond this region, the charge distribution is dominated 
by screening and falls off in a DH-like manner after approximately one Debye 
length, as it is shown in the inset of the figure [TOl 

At low density (n < O.OIM), ri„ depends on ^ in a similar way to what is shown 
in figure [151 For n = O.OOIM, we can associate the smallest value of ri„ to 
the Manning radius within the range 1 ^ C ^ 3. Besides, that range roughly 
discriminates the region where ^ren passes from an increasing to a decreasing 
function, as it can be seen in figure [T71 When ^ > 3.4, infiection points related 
to the finite size of microions enter to scene, indicating the accumulation of a 
second layer of counterions and consequently, the additional decrease in ^ren- 
Further increasing ^ produces annihilation of the inflexion points at ^ w 4 and 
simultaneously go+{r) starts exhibiting a peak, as flgure [TOl illustrates. 

Previously, when the charge inversion phenomenon was discussed, we associated 
the onset of a peak in the coions density profile to the settling on of a "layered" 
ionic cloud and we argued that the renormalized charge is determined by the 
net charge enclosed in that cloud (see figure [TB)) . A similar conclusion has been 
given from the comparative study at very low ionic strength of the electrokinetic 
and the effective charges, which were obtained from electrophoresis and light 
scattering data of suspensions of latex particles, respectively pT| . 

According to the study carried out in it may be stated that, interacting 
at large distance, colloidal particles act as if they carry an effective charge 
equal to the one enclosed by the so-called outer Helmholtz plane (OHP), from 
which the diffuse double layer begins. To a certain extent, the definition of 
OHP can be associated to the choice of the r value, from which go+if) and 
go-{r) fall off symmetrically with respect to 1 in a DH-like form, i.e. the choice 
of an adequate range to fit correlation functions within the linearized theory. 
On the other hand, we would like to point out that the determination of the 
electrokinetic charge depends on the previous evaluation of the zeta potential 
and, as it was aforementioned, the behavior of ^^en shown in figure [T7] and that 
of the zeta potential for a charged cylinder found in |18] are qualitatively similar. 
In this way. our results suggest that the HNC approach captures appropriately, 
at least qualitatively, the short-range correlations determining the double layer 
structure. 
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On the other hand, when the macroion is immersed in an asymmetric electrolyte 
composed by divalent and monovalent ions, the renormalized charge evaluated 
from HNC calculation behaves as shown in figure [201 Again, the overshooting 
and charge inversion phenomena are observed. In the 2:1 case at n = 5 x 10~^M 
(k(Tc ~ 0.09) and n — O.OOIM (kctc « 0.13), figure suggests the presence of 
a well defined layer of condensed counterions. At these densities, the Manning 
radius dependence on ^ is qualitatively similar to the one predicted by PB 
theory; however, it appears to underestimate the thickness of the condensed 
layer. If we would consider the renormalized charge as the net charge enclosed 
by the Manning radius (two-state model), then the renormalized charge should 
be smaller than the bare one, and therefore, the overshooting phenomenon would 
not be present. 

In spite of a thicker condensed layer, depending mainly on short-range cor- 
relations between the macroion and the counterions, the HNC approach also 
predicts overshooting. This phenomenon lies in the fact that divalent colons 
are expelled further away in the double layer than the monovalent ones in the 
1:1 situation, which results in a stronger electrostatic potential at larger dis- 
tance [H] . Although the PB theory accounts for this phenomenon, figure [20] 
shows that the mean field solution underestimates the magnitude of the over- 
shooting effect as well as the range of ^ where it is present, which is particularly 
notorious at high density. 

Likewise, at high enough density (for instance, n = 0.005M (kCTc w 0.29)), as ^ 
increases the Manning radius increases faster that its PB counterpart (dashed 
line in figure [2T|) and eventually fuses with the inflection point characterizing the 
salt. Again, when inflection points are close enough, a peak in go+{r) appears, 
i.e. a local increase of colons density takes place, as it can be seen in the bottom 
frames of figure 1221 In these circumstances, the Manning inflection point and 
the Debye length become comparable so that, the charge renormalization is 
determined by a more correlated ionic cloud giving place to a saturation value 
of Cren within the studied range of ^, as it is displayed in figure [201 

With respect to the 1:2 case, the differences between PB and HNC results 
become more significant. For all the considered densities, soon after the linear 
charge density is larger (in absolute value) than the Manning condensation 
threshold {^m — —0.5), the renormalized charge decays quickly to zero and 
gets the positive sign at ^ < —1.3 (see the inset of figure 1^0)) . As well, there 
is no qualitative similarity in the Manning radius evaluated from both of those 
approaches, suggesting that the short-range correlations between macroions and 
microions play a more predominant role in this situation. Although in figure 
the colon densities for both 1:2 and 2:1 electrolytes exhibit a qualitatively similar 
behavior as |^| increases, we call the attention on the difference between their 
orders of magnitude at r ~ cjc- in the 2:1 case the first microion layer around the 
macroion corresponds only to counterions whereas in the 1:2 case, this appears 
to be conformed both of counter- and colons. This fact indicates a change in 
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the short range behavior of the mean force potential between the macroion and 
the coions, as it was mentioned earher. In the 1:2 case, increasing further |^|, 
the counterion fraction of the ionic cloud around the macroion does not only 
compensates for its total charge but even exceeds it, giving rise to a charge 
inversion phenomenon, as it is displayed in the inset of figure 1201 

Also, the results show the existence of a maximum value of the effective charge 
In order to describe the dependence of that maximum on the macroion 
radius (Tq, we evaluate ^^en at different values of (Tq for both symmetric and 
asymmetric electrolytes. The outcome of these calculations within the HNC 
convergence range are shown in figures [23l and [24l As it can be seen, the qual- 
itative behavior appears to be the same previously discussed: the observation 
of the maximum of S,reniO followed by a monotonic decrease is systematically 
found in all studied cases. In figure 1251 the dependence of the maximum of 
^ren On thc rcduccd radius kCTc [o'c = o"o + is displayed. Although system- 
atically lesser, £,™n^ qualitatively follows the same tendency (at less within the 
considered range) that the saturation value from PB theory (^f^^) at kCTc > 1 
in both symmetric and asymmetric electrolyte [ini [H] . We can associate the 
Crln^ to the counterion condensation in the sense that, around the bare charge 
^* at which thc maximum is reached, the renormalized charge is mainly fixed 
by the macroion-counterion correlation. In this sense, since we use a 2D system, 
which overestimates correlations between microions, figures [25] can be used to 
establish both upper (PB theory, no microion correlations) and lower bound to 
the exact value of . 

Finally, coming back to the figure [23] for n = O.OIM and ctq — 25A, a particular 
behavior of S^ren is found as ^ > 15. At such conditions the decrease of £,ren 
appears to be less pronounced that the one found for 8 < ^ < 15. This fact 
would indicate the eventual appearance of a plateau region in ^ren at large 
enough ^; nevertheless, since the HNC calculation does not converge at Z > 16, 
it is difficult to give a conclusive argument in this direction. 



5 Conclusions 



We have evaluated the renormalized charge in a 2D system formed by a central 
charged particle (macroion) immersed in an electrolytic solution. By solving 
numerically the Ornstein-Zernike equations within the framework of the hyper- 
netted chain approximation (HNC), we have calculated the correlation functions 
between the macroion and the microions {go±(j'))- Under the assumption that 
the correlation functions fall off exponentially, we found the renormalized charge 
by fitting the asymptotic behavior of the correlation functions to the one de- 
rived from the linearized PB theory (Debye-Hiickel limit) . The adequacy of this 
procedure was proved by comparing the numerical results to the analytical ones 



16 



from exactly solved 2D models. 

The 2D model could be applied to study a cylindrical macroion immersed in an 
electrolyte formed by cylindrical co- and counterions, for instance, charged stiff 
polymers. Moreover, since particles in the model interact through a logarithmic 
potential, it would be expected that it preserves aspects concerning the nature of 
the Coulomb interaction, i.e., in a 3D model, the dependence of the renormalized 
charge on the bare one would share the same qualitative characteristics that 
those described here. 

Recently, by using the theory of integrable systems, some analytical results con- 
cerning the renormalized charge within the non-linearized Poisson-Boltzmann 
(PB) theory have been reported [TU]. The results of this work show that the 
renormalized charge is not explained qualitatively by the PB theory. As the 
bare charge of the macroion increases sufficiently, the PB theory predicts a con- 
stant value of the renormalized charge; on the contrary, the HNC renormalized 
charge reaches a maximal value and then monotonically decreases, both in sym- 
metric and asymmetric electrolytes. This behavior takes place because the HNC 
approach appreciates the difference between the potential and the potential of 
mean force, which is exactly what the Poisson-Boltzmann theory neglects. 

On the other hand, the structure of the ionic cloud around the macroion, which 
provokes the charge renormalization, can be described through either the den- 
sity profiles go±{f) or the inflection points of the integrated charge density P(r) 
(defined by eq. (fT7|) ). For a weakly charged cylinder, the short range behavior 
of the potential of mean force between macroion and microions appears to corre- 
spond to the Coulomb interaction. Also, there is one inflection point associated 
to the screening phenomena in the bulk of the salt. Consequently, at that limit, 
the density profiles suggest a diffuse cloud, which gives rise to a linear screening 
that can be described by the Debye-Hiickel approximation. 

At low density, as the bare charge ^ increases, the stronger accumulation of 
counterions around the central particle results in a change of the short-distance 
behavior of the potential of mean force between macroion and colons, which has 
been interpreted as a precursor of the (Manning) counterion condensation [5]. 
Simultaneously, according to the inflection point criterion, a new inflection point 
at rm < can be identified. This point can be associated to the Manning 
radius provided that the bare charge is close to the Manning threshold; in such 
situation the renormalized charge is mainly fixed by the macroion-countcrion 
correlation. Besides, since the maximum of the renormalized charge (C^n^) 
exhibits the same characteristics as the plateau of the PB theory, it can be 
assumed that under these circumstances, the thermal energy of a counterion 
is nearly compensated by the reversible work required to remove it from the 
vicinity of the macroion. 

On the other hand, by studying $,"1^^ within an intermediate range of reduced 
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radius (10~^ < kctc < 1), we can conclude that ^"^^"^ as a function of kCTc is 
systematically smaller and present the same qualitative behavior than the satu- 
ration value of the PB theory. Since in our model the microions interact through 
a logarithmic potential, the role of microion correlations has been overestimated 
from the point of view of 3D systems. In this way, our HNC results could be 
established as lower bounds for the real ^l^n^. 

As ^ increases further, r„i becomes bigger until it fuses with the inflection point 
characterizing the electrolyte. Under such circumstances, most of the counteri- 
ons are located at the cylinder surface, whereas the rest of the screening charge, 
is spread in the ionic cloud. Nevertheless, the merging of the inflection point is 
related to the onset of a peak in the colons density profile. This suggests that 
beyond the countcrion layer, the ionic cloud is not a diffuse layer which can be 
handled by means of the linearized theory, as the two-state model claims, but 
a more complex structure is settled by correlations between microions. 

If we consider an asymmetric 1:2 electrolyte, the net charge in the vicinity of the 
cylinder surface has opposite sign to that of the bare charge. This fact causes 
that the colons at some distance from the surface to be attracted towards the 
macroion (giving rise to the mentioned peak) and that the counterions to be 
repelled from the macroion, leading to a reversion of the charge distribution in a 
region some distance from the cylinder. In this way, the evaluation of the renor- 
malized charge from the asymptotic behavior of go±(r) allows to account for 
the charge inversion phenomenon without additional assumptions, in contrast 
to PB theory, in which the introduction of the countcrion binding conjecture 
becomes necessary. 

Therefore, the present study suggest that the HNC approach is able to take 
into account, at least qualitatively, the behavior of correlations responsible for 
the countcrion condensation process, as well as the charge inversion one. As 
well, it gives indications in order to obtain a better description of the charge 
renormalization, because it intrinsically explains for the structure of the electric 
double layer. So, a theoretical method to estimate more adequately the renor- 
malized charge at high enough ^, must consider a more detailed description of 
the electric double layer structure by taking into account short range correla- 
tions between macro- and microions and between the last ones (for instance, 
density functional or dressed-ion theories). 

Finally, a question remain to be solved: Does ^ren always continue to decrease 
at large ^? A unique of the outcomes presented in this work appears to indicate 
the possibility of a plateau region in ^ren at large enough f ; however, due to 
the non convergence of HNC calculation at high ^, wc cannot categorically state 
this conclusion. As it was mentioned, HNC approach is inadequate to treat the 
ion clustering effects, which could be of extreme importance at such a strong 
coupling. A device to overcome this difficult might be the introduction of a 
bridge function, which takes indirectly into account correlations mediated by 
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two or more particles. 
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Figure 1: Correlation functions go±{r) at different values of ZT for 
a guest charge immersed in a symmetric TCP with coupling para- 
meters r = 0.2 (top) and T = 0.4 (bottom). 
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Figure 2: Correlation functions between microions of the symmetric 
TCP at different values of coupling parameter T. The plotted func- 
tions correspond to the ones evaluated when ZT — 3; however these 
functions are independent of the magnitude of ZT, as expected from 
the equations (13])-(15]). 
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Figure 3: Functions Lo±(7') used to evaluate renormalized parame- 
ters. The curves were calculated from the correlation functions at 
ZT = 3. Triangles delimit the r range used for the linear fit. 
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Figure 4: Renormalized charge as a function of the bare one for 
different couphng parameters. Dashed and thick continuous hues 
represent the mean-field PB analytical result for cto = ct = and the 
numerical one when kctq = 10^^, (7 = 0, respectively [5]. Dotted line 
indicates the exact result from field theoretical arguments (kctq 
0,(7 = 0, r = 0.5) U. 
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Figure 5: aos factors characterizing the short-distance potential of 
mean force at T = 0.1 and F = 0.5. In order to exclude the hard-disk 
effect, the data shown was calculated from a linear fit of In g(r) at 
6cr < r < lOcr « O.Ik^^. Dotted lines represent theoretical limits 
(HH) at F = 0.5. 
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Figure 6: Integrated charge distribution as a function of the radial 
distance at F 0.4 and different values of ZT. Dashed lines repre- 
sent the distribution of counterions when ZT = 0.4 and ZT = 3.6. 
Radial range corresponds to two Dcbyc lengths. 
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Figure 7: Inflection points of tlie cfiarge distribution as a function of 
zr for different values of F (left). Continuous line correspond to the 
numerical solution of PB equation (data taken from [TO])- Detailed 
plot showing the Manning radius (right). In the inset are shown 
the correlation functions go-{r) at F = 0.1 for ZT = 2,3,4, which 
suggest that the onset of a well-defined second layer of counterions 
provoking smaller values of ri„. Thin dashed lines indicate the po- 
sition of one, two and three particle diameters. Dotted lines are a 
guide to the eye. 




Figure 8: Correlation functions between microions at F = 0.2 (left) 
and F = 0.7 (right) for the asymmetric TCP. Dotted lines denote 
theoretical asymptotic limits given by field theoretical arguments 
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Figure 9: Correlation functions go±{r) for a guest charge immersed 
in an asymmetric TCP for different values of Z > (1:^ case). The 
curves were evaluated at F = 0.2 (top) and F = 0.7 (bottom). 
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Figure 11: Functions Lo±{r) used to evaluate renormalized pa- 
rameters. The curves were calculated from the correlation func- 
tions at r = 0.7. Triangles delimited the range of r according to 
10-^ < I log^o. Wl < 10"^ for- ^ ^ -4. 
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Figure 12: Renormalized charge for a guest charge immersed in an 
asymmetric TCP at different couphng parameters. Dotted line in- 
dicates the condition Zren = Z. In the inset the HNC result for 
r = 0.9 is compared to the one predicted by PB theory (dashed 
line). 



30 




r/2<j r/2a 



Figure 13: Short distance behavior of correlation functions go±ir) 
for Z = —3.6 and different couphng parameters. 
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Figure 14: aos factors characterizing the short-distance behavior of 
the potential of mean force for different values of F. Dotted lines in- 
dicate factors related to the Coulomb potential (|ao^"'^| = |Zz±r|). 
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Figure 15: Inflection points of charge distribution as a function of Z 
and r for ^il (left) and 1:^ (right) cases. At the left, the inset shows 
an additional inflection point at larger distance when T = 0.8, which 
is related to the charge inversion. Dotted lines were obtained by 
interpolation but are intended as a guide to eye. As in the symmetric 
case, the smallest values of are related to the finite radius of 
particles, as dotted lines indicate in the inset at the right. 
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Figure 16: Absolute value of the charge density around guest particle 
for the ^:1 case at F = 0.8. Peaks indicate a change of sign of p{r) 
and therefore a change of sign in the renormalized charge. Note that 
the onset of such peaks imposes a restriction on the r-range used to 
the fitting procedure. 
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Figure 17: Dependence on the salinity of the renormahzed charge 
for a charged cyhnder immersed in a 1:1 electrolyte {[n] ~ M). The 
renormalized linear charge density is compared to the exact PB re- 
sults for n = O.OOIM (kctc ~ 0.074, dotted line) and n = O.IM 
(kcTc ~ 0.74, dashed line). 
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Figure 18: Inflection points in the integrated charge distribution 
(left). Arrows indicate values of kctc for each density {[n] = M). 
Dotted lines are intended as a guide to the eye; the continuous line 
indicates the exact PB resuh for n = 0.001 {kcfc ~ 0.074) [lO]. At 
the right are shown the correlation functions go± for ^ = 3 — > 6 at 
n = 0.001 M. At higher values of ^, the correlation functions exhibit 
inflection points as a consequence of finite radius of microions, which 
suggest the formation of a second layer of counterions. 
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Figure 19: Correlation functions 5o±('') at n = O.OOIM (top) and 
n = O.IM (bottom) for a cylindrical macroion. Dotted lines delimit 
the region r = CTc + 2cr and r = CTc + 4(t. When n = O.OOIM and 
^ = 4, the onset of a peak in go+(^) takes place. This point is 
associated with the annihilation of inflections points (see the figure 
(fT8|) ). In the inset is shown the charge density at n = O.IM for 
^ = {1, 2, 3, 4, 5} (dashed line indicates ^ = 3). 
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Figure 20: Rcnormalizcd charge for a cylindrical macroion im- 
mersed in an asymmetric electrolyte at different salinities. Dotted 
and dashed lines indicate analytical PB solutions at ti = 0.0005M 
(kCTc ~ 0.9) and n = O.OIM (kCTc ~ 0.41), respectively. 
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Figure 21: Locus of the inflection points of P{r) for a charged cylin- 
der of radius (Jq = lOA immersed in an 1:2 (left) and 2:1 (right) 
electrolyte. Continuous and dashed lines correspond to calculation 
of by solving the PB equation for nac = 0.1 and nac = 0.3, 
respectively [TO]. Note that for n = 0.0005M and n = 0.005M corre- 
spond KCTc ~ 0.09 and kctc ~ 0.29, respectively. Arrows indicate the 
K(Tc limits. 
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Figure 22: Correlation functions gQ±{r) for 1:2 (top) and 2:1 case 
(bottom) at n = 0.005M. In the inset are shown the associated 
charge densities (dashed hne is used for ^ = —1.0). Dotted lines 
marks r = CTc + 2(t and r ~ ac + 4cr. 
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Figure 23: Rcnormalizcd charge in the symmetric case for different 
values of colloid radius ([cto] = A) at n = O.OOIM (left) and n ~ 
O.OIM (right). Dotted lines indicate the PB solutions from [16] for 
ctq = 150a (kctc ~ 1-6 and 3.6, respectively). Arrow indicates a value 
of ^ from which a less pronounced decrease is apparently found. 
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Figure 25: Maximum value of renormalized charge as a function of 
the reduced radius kctc for symmetric (left) and asymmetric (right) 
electrolyte. At the left, continuous and dashed lines arc the theo- 
retical PB upper bound for Zj-en at saturation (analytical prediction 
from [TOj accurate for large a) and the PB saturated charge evaluated 
numerically in [10] . respectively. At the right, dashed and contin- 
uous lines are, respectively, the PB saturated charges for 2:1 and 
1:2 cases [12]. Arrows indicates theoretical PB limit when Ka^ 0. 
Dotted lines are intended as a guide to eye. 
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